JOURNAL OF RESEARCH of the National Bureau of Standards— B. Mathematics and Mathematical Physics 

Vol. 67B, No. 2, April-June 1963 

Effect of Molecular Weight on Viscoelastic Properties of 
Polymers as Predicted by a Molecular Theory 

Hansjorg Oser and Robert S. Marvin 

(February 27, 1963) 



Calculations have been made covering the predictions of a model representing the vis- 
coelastic behavior of rubberlike polymers for molecular weights greater than Mc, Mc being 
the lower limit of the range in which the viscosity is proportional to M^'^. A pronounced 
difference in the character of G" is predicted for polymers with molecular weights b(^tween 
five and ten times Mc as compared with those whose molecular weights arc outside this range. 



In a recent publication [1] ^ we presented the steady 
state and transient behavior predicted by a model 
designed to represent the mechanical behavior, in 
shear, of rubberlike polymers. This model is very 
similar to the molecular theory of Rouse [2], except 
that it includes du'ectly the effect of entanglements 
found in long polymer chains. 

Our model assumes an entropy-type elasticity as 
the origin of the relaxing force observed in a stress- 
relaxation experiment. It involves a resistance 
coefficient per unit segment of a polymer chain 
which is evaluated from the steady state viscous be- 
havior, and the molecular weight of the polymer, 
expressed as the ratio of molecular weight to critical 
molecular weight, M^ defined as the lower limit 
of the range in which the viscositv is proportional 
toMK 

For the particular example we presented in 
reference 1 tliis ratio was 60, and at the time we 
prepared this previous paper wt did not have system- 
atic calculations showing the predictions of the 
model for a range of molecular weights. We made 
an effort (sec. 3, p. 176 of reference (1)) to relate the 
extrema of G^^ and J'^ to the molecular weight 
through the use of various approximate relationships, 
since the full expressions can only be treated numeri- 
cally. We have now carried out further numerical 
calculations covering a range of molecular weights, 
and jfind that our previous approximate relations 
were completely inadequate to represent the posi- 
tions of these extrema as functions of molecular 
weight. In addition, we have found certain quite 
unexpected qualitative features of the predicted 
curves when the molecular weight is varied between 
one and ten times the critical molecular weight. 
Since these qualitative features may present the 
possibility of a rather critical test of the validity 
of this model, we present the results here in some 
detail. 

The results of our previous paper may be rep- 
resented in terms of normalized functions, the 



normalized modulus, 

and the normalized compliance, 

3*(7)=7J*=3'W-i3"W. 

C2 is a parameter equal to the kinetic theory equi- 
librium compliance of a vulcanized rubber with 
molecular weight Mcl2 between cross links. Obvi- 
ously, these are reduced functions since the ordinary 
temperature and density dependence are included in 
C2. T and 6 are reduced frequencies, defined by 



^=B{l-\-i)—'yJiu)C2r 






(1+i) 



where co is the actual radian frequency and r the 
viscosity of a polymer of molecular weight Mc. In 
these terms, 



@*(7) = 



tanh 7+L tanh Ky _ 1 
l+L tanh 7 tanh Ky~^ 



(1) 



where 



K=Lil-l). 



l\l/2 



M 



1 Figures in brackets indicate the literature references at the end of this paper. 



Equation (1) represents the predictions of our model 
up to frequencies where the modulus is within one or 
two decades of the limiting high frequency modulus. 
This range includes all the extrenui which depend on 
molecular weight. The maximum in (?'' occurring 
just before G' reaches its limiting or glassy value is 
independent of molecular weight. Its position and 
magnitude can be calculated easily as given in the 
previous paper [1]. 

The calculated values of the normalized modulus 
and compliance as a function of the normalized 
frequency, 6, are shown for several values of Z, used 
here as a measure of molecular weight, in figures 
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Figure 1. Real part of dynamic modulus for various molecular 

weights. 

Parameters are values of log Z=log (M/Mc). 
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Figure 2. Imaginary part of dyriamic modulus for various 
molecular weights. 

Parameters are values of log ?=log (M/Mc). 
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Figure 3. Real part of dynamic compliance for various mo- 
lecular weights. 
Parameters are values of log Z=log (M/Mc). 
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Figure 4. Imaginary part of dynamic compliance for various 

molecular weights. 

Parameters are values of log Z=log (MjMc). 
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Figure 5. Reduced frequencies of the extrema as a function of 

1 = ^^ 

Me 



1 to 4. We see that ®'' and 3'' do not show an}' 
extrema for certain values of I, and that for other 
values of I there are two minima and two maxima 
in ®" below the glassy range. 

This occurrence of two sets of extrema was cer- 
tainly unexpected. To illustrate this behavior more 
clearly, we show in figure 5 a plot of ^ext, the fre- 
quencies at which the extrema occur, as a function 
of I and L. We distinguish four separate curves 
A and C belonging to ^" (6) , B and D belonging to 
®'' {6) . The lower branches of A and C represent 
at the minimum of 3'', the upper branches 6 at the 
maximum of that function. The lower branches of 
the curves B and D represent 6 at the maximum of 
®" and the upper branches give d at the minimum 
of that function. 
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The points with vertical tangents (fig. 5) corre- 
spond to inflection points of ®'' or 3'' with a horizontal 
tangent. This follows immediately, once we realize 
that they represent points where a maximum and 
mininmm have coalesced. 

It is interesting to examine the asymptotic behavior 
of the branches of C and D^ for l-^ ^ , since here these 
curves apparently become straight lines and we ex- 
pect simple power laws for the functional relation 
between Q and I (or L). We observe from figure 5 
that a common lower bound for all four branches of 
the curves C and D is given by 



1 



21' 



1 



as indicated by the dotted line in figure 5. Hence 

in the region /> 10, Z>16.5. 

This inequality shows that tanh {Ky) in the ex- 
pressions for (y and 3^''can be put equal to unity and 
therefore does not influence the location of any of the 
extrema for large /. Our functions simplify then as 



follows, if we introduce x = tanli 



G) 



@''(6') = Im{7COth (7+x)}=^ 



3-(^)^_In,/itanh(7+x)\ 



sinh 2((9+x)— sin 2(9 
cosh 2(i9+x)— cos2i9 

(2) 



1 sinh 2((9+x)— sin 2(9 
'2i9cosh2((9+x)-fcos2^ 



(3) 



We turn now to the determination of the 
asymptotic behavior of the branches I, II, III, IV. 
On both branches I and II of curve D, 6 and x approach 
zero as L increases. We can thus use Taylor's 
expansion to obtain: 



&''=-. 



ex 



1 



e{2e+x)+0{e^). 



2d^+2xe+x' ' 3 

The derivative of the first term (with respect to 6) pre- 
dicts a maximum at 6^^^=-^^=-^^ ( X ~y for large L j 

with ©max =0.207. This is the only extremum pre- 
dicted by the first term alone. To find a condition 
for the minimum we re-write the above as: 

this expression being permitted because on branch 
II X^d^. This latter relation justifies our taking, 
for large L, 



©' 






From which we predict a miniinuni at 



with a magnitude 






An expression for the asymptote of the minimum 
in 3^' represented by branch III (curve C) may be 
obtained from a similar expansion of tanh (7+X), 
again retaining two terms and using the fact that 
X<if. It yields, as is apparent from figure 5, the 
same value as the asyniptote for branch II: 



from which 



and 



^ 26/^3 ' 






x:n=i.04z- 



The maximum in 3''? branch IV of curve C, 
appears to be nearly independent of L If it were 
strictly independent, it would be the maximum of: 

J_ sinh 2^— sin 26 
2^*cosh26>+cos26>' 

the asymptotic expression corresponding to eq (3) 
for large /. This function has a maximum value of 



at 



Xax-0.4172 



dm 



:1.127. 



The following tabulation shows that these maxima 
of our calculated 3" functions depend only very 
slightly on I: 

for Z=5.012: X'a.=0.4325, ^„,ax=1.019 

for Z-100.00: 3:;.= 0.4178, ^„,,-1.122. 

It is difficult to do more tlian point out the exist- 
ence of curves A and B (fig. 5), which was not even 
suspected prior to these detailed numerical calcula- 
tions. Since in this region tanh Ky cannot be taken 
as unity, the analytic expressions for ©'' and 3'' 
are very complicated. There is no reason to expect 
that expressions for curves A and B would be simple 
enough to be useful. Therefore we leave them with 
the presentation of the numerical results. These 
computations have now been carried out in sufficient 
detail that we are essentially certain we have found 
all the extrema predicted by the model for both 
®''and3''. 
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The behavior of ®' and ©'' calculated here as a 
function of molecular weight should provide a 
critical test of our model if appropriate experimental 
data could be found for comparison. In the molecu- 
lar weight range of about 5 to 10 times Mc this model 
predicts an essentially constant value of G^ ' over two 
to three decades of frequency, with Q' increasing 
fairly steadily over the same range. For higher 
molecular weights the G" curve should show a pro- 
nounced maximum and minimum, their separation 
becoming greater on both the frequency and modulus 
scales with increasing molecular weight. In the 
vicinity of the minimum of G'\ G' should remain 
fairly constant. Since we do not know how a 
distribution of molecular weights should be intro- 
duced into our model, it would be preferable to 
make the comparison with experimental measure- 
ments on fractionated materials, but lacking such 
data it seems most reasonable to use a weight 
average molecular weight. 

While we do not know of any dynamic measure- 
ments which can be used to check these predictions 
directly, the stress relaxation measurements and 
calculated relaxation spectra reported by Tobolsky 
and Murakami [3] for a series of polystyrene samples 
with weight average molecular weights ranging from 
4.25X10^(-1.1 Me) to 3.15X10'(=8.3 M,) (theh 



sample numbers 6 to 15) and MJMn=l^5 represent 
a reasonable check on the predictions summarized 
above. The stress relaxation modulus is approxi- 
mately the mirror image of G\ and the spectrum 
H{t) is, to a zero order approximation, equal to 
(2MG"il/r). _ 

The graphs of H shown in figure 4 of reference 3 
show no plateau region for molecular weights between 
Mc and 2.4 Mc (samples 6 to 11). Distinct plateau 
regions are seen for samples with molecular weights 
between 3.2 Mc and 8.3 M^ (samples 12 to 15). For 
all these samples the G{t) curves are monotonically 
decreasing, with no hint of a plateau region. For 
higher molecular weights, of course, we normally find 
a region where G(t) or G^(o:) are changing so slowly 
as to suggest strongly a plateau. This type of be- 
havior is observed in our calculated curves for Z>15. 
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